Optimizing inhomogeneous spin ensembles for quantum memory 
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We propose a method to maximize the fidelity of quantum memory implemented by a spectrally 
inhomogeneous spin ensemble. The method is based on preselecting the optimal spectral portion of 
the ensemble by judiciously designed pulses. This leads to significant improvement of the transfer 
and storage of quantum information encoded in the microwave or optical field. 

PACS numbers: 03.67.Lx, 42.50.Pq, 76.30.Mi, 85.25.Hv 



I. INTRODUCTION 



Recent experimental demonstrations of strong cou- 
pling between spin ensembles (SEs) and microwave pho- 
tons of superconducting resonators [lUff are an impor- 
tant step towards realizing functional, hybrid quantum 
devices [7H18|. Such hybrid devices may benefit from 
combining the advantageous properties of very different 
subsystems, or "blocks": (i) a quantum processor block 
containing e.g., q uantum dots [19|, [20] or superconduct- 
ing qubits |2l| - |24| which can perform rapid quantum gate 
operations but are vulnerable to decoherence due to their 
strong coupling to the environment and/or the noise of 
the external controls; (ii) a quantum memory block con- 
sisting of a SE of active dopants in a solid [lUs [25|, [26| or 
trapped ultracold atoms J7H15| which are weakly coupled 
to the environment and therefore are suitable for infor- 
mation storage; (hi) a quantum "bus" or interface, such 
as a microwave cavity, whose interaction with the other 
blocks can be quickly switched on and off by, e.g., tun- 
ing in and out of resonance |27H3fj| . A related scenario 
concerns reversible transfer and storage of optical exci- 
tations in SEs used as memories for photonic quantum 
repeaters [3lM37| . 



The potential advantages of hybrid quantum devices 
are countered by decoherence during the transfer and 
storage of quantum information (QI), which is rooted 
in the homogeneous (lifetime) broadening and the in- 
homogeneous spectral width of the SE constituting the 
memory block. Using magnetic dipole or optical Raman 
transitions can greatly prolong the relaxation time of the 
spins. But since typically the coupling of individual spins 
with the microwave cavity or optical field is weak, QI 
must be collectively encoded in an ensemble of N 3> 1 
spins, which are prone to spectral broadening J18l 13a. [39| . 
This broadening introduces dephasing on the memory 
block which can both hinder the transfer of QI and limit 
its coherent storage time. Although conventional spin- 
echo techniques [40|-|43| can prolong the storage time, 
they may be impractical for spectrally broad SEs, and 
they cannot improve the QI transfer (in fact they would 
preclude the transfer altogether). 



Here, after analyzing the fidelity of SE quantum mem- 
ory, we present a method that can significantly improve 
its performance. The method is based on preselecting the 
optimal spectral portion of the inhomogeneously broad- 
ened ensemble using judiciously designed pulsc(s), fol- 
lowed by the transfer of all the spurious spins to an aux- 
iliary metastable state. Although the selected subensem- 
ble contains far fewer spins and therefore requires longer 
QI transfer times, we find that the resulting fidelity of 
the quantum memory during the QI transfer and storage 
is much higher under experimentally realistic conditions. 



II. FIDELITY OF AN ENSEMBLE QUANTUM 
MEMORY 



To set the stage for the discussion, assume that a 
qubit is encoded in a SE as a superposition of its col- 
lective ground state \ipo) — \gi,g2,.-.,gN) an d the 
fully-symmetrized single-excitation (Dicke) state \ipi) = 
N ~ 1/2 Y,j\J), wh ere \j) = \gi, g 2 , ■ ■ ■ ,e h . . . , g N ) de- 
notes a state with only spin j excited. If all the spins 
had the same resonant frequency loq on the transition 
\g) — > |e), the single-excitation state would evolve in 
time r > as \4>i{t)) = iV" 1 / 2 ^- e~ iLU " T \j), while the 
ground state l^o) remains unchanged. Yet, due to in- 
homogeneous broadening, each excited spin has differ- 
ent resonant frequency ujj . Then, even if the symmetric 
state |t/>i) is prepared at r = 0, it would evolve into 
\4'i(t)) = 7V-V2 y\ e -«^-r |i). We may thus define the 
storage fidelity at time t > as the squared overlap of 
state | - 0i(r)) with its inhomogeneously-broadened coun- 
terpart |i/'i(t)): 



F(t)^\(Mt)\Mt))\ 2 = 



-^ /"n(w)e~ 4( ^ o)T du; 



(1) 

where n(w) is the ensemble spectral density normalized 
to the total number of spins, / n(oj)du = N. Note that 
-F(O) = 1 for any n(w). Since the ground state IV'o) does 
not evolve in time, F(t) quantifies how the fidelity of in- 
formation already encoded in a SE decreases over time. 



Specifically, for an ensemble with Lorentzian spectrum 
of width A, n(oj) = n [l + (u - Wo) 2 /^ 2 ] -1 , the fidelity 
loss is exponential in time, F(t) = e~ 2Ar . As shown be- 
low, Eq. ([1]) also characterizes the efficiency of excitation 
transfer to and from the SE. 



A. Transfer fidelity 

Consider a single-mode field of a microwave cavity 
or an optical beam interacting with the SE. Assume 
that initially the field contains a single excitation (pho- 
ton) of frequency wo and the SE is in the ground state 
IV'o}- Each spin j in the ensemble interacts with the 
field on the transition \gj) — > \ej) with the coupling 
strength r\j . In the rotating wave approximation, neglect- 
ing the field and spin relaxations, their combined state 
at any time r can be written as |\&(t)) = a(r) \l,ipo) + 
EjftW^'^^'lO.i)- where |1, -0o) refers to the 
state with a single photon and all the spins in the ground 
state, and |0, j) denotes the field vacuum and the jth spin 
excited. The probability amplitudes a and j3j evolve in 
time according to 



a 



: E^ft 



-i(ojj —ojq)t 



fa = -irisae'to-" *, 



which yields d(r) = — V. |?7j| 2 / dr'a(r')e" 
Assuming that the spin-field couplings 77^ are not corre- 
lated with the transition frequencies uij [3j-l5(, we finally 
obtain 



d(r) 



-fj 2 N 



dr'a(r')V^( T - T ')> 



(2) 



where rj 2 = N 1 J2 t \rji\ 2 and F(t) is given by Eq. ((T|). 
With F(t) = 1 for all r, Eq. © predicts Rabi os- 
cillations between the field and the SE according to 
a{r) — cos^v/Vt), so that at time r tr = ir/fjyN there is 
a full retrieval of the excitation into the field, |a(rt r )| 2 = 1 
[with the combined state |\&(Tt r )) = — 1 1 , V^o ) acquiring a 
■k phase shift] . In the presence of inhomogeneous broad- 
ening, however, F(t) decreases with time, resulting in 
damped Rabi oscillations. The fidelity of the transfer 
followed by retrieval is the value of |a(rt r )| 2 after one 
such Rabi oscillation. For F(t) < 1 the fidelity loss is 
0(1 — F). Hence, increasing the storage fidelity F(t) for 
all times r G [0,7t r ] also improves the transfer fidelity. 



B. The optimal spectrum 



of the subenscmblc. To this end, we may consider two 
different tasks: (a) QI transfer to and from the memory, 
and (6) QI storage in the memory for a specific time r s . 
(a) For QI transfer, we require that the fidelity F'(t) 
for the selected subensemble be high for all r G [0,r tr ]. 
For any symmetric spectrum n'(uS), Eq. (JlJ can be ex- 
panded in a Taylor series 



F'(r) = 



E 

fc=0 



(2fc)! 



k^2k 



-<(w-wor> 



l-iiu-LJofy, 



(3) 

where ((■■•)) = / n'(cu)(- ■ ■ )duj/ J n'(uj)duj denotes the 
average over the spectral distribution n'{u>). Hence, the 
spectral variance ((w — wo) 2 ), which is the leading term 
of the expansion, should be as small as possible for a 
required number N' of the selected spins. We then find 
that the optimal spectrum is n'(uS) = n(u) for \u — luq\ < 
A and n'(cu) = otherwise, where A is such that N' = 
Iuj°-a n i UJ )doJ. In other words, we should select all the 
spins from the frequency interval ui G [wo — A, wo+A], and 
none outside of it. Assuming an approximately constant 
original spectrum n(w) ~ n around to = cj , we obtain 
N' ~ 2An and 



F'{t) -sine 2 (At) ~ I 



|A 



2 T 2 



(Ar<l). (4) 



-i(uij— « )(t— t') 



(6) For QI storage over time r s , we require that the 
fidelity F'(t) be high at exactly t = t s . For a symmetric 
spectrum n'(w), the fidelity of Eq. ([I} reads 



F'(r s ) = 



1 

A 7 



n (w) cos([w — u)o]T s )duj 



(5) 



Hence, we should select all the spins with frequencies uj 
which maximize cos([u; — ujq\t s ), i.e., a comb-like spec- 
trum n'(uj) peaked at ui ~ u>k = luq + 27rfc/r s , where 
k = 0,±1,±2, ... For a required number of selected 
spins TV', the optimal spectrum is then composed of a 
series of rectangular distributions around frequencies uik- 
n'{ui) = n(w) for \u — u>k\ < A and n'(ui) = otherwise 
(we note a similar result in |33|). Assuming the original 
spectrum changes little around each u = Wfe, n(oj) ~ rik, 
we obtain N' ~ 2A J^k n k an d 



F'(t b ) 



smc 



! (Ar s ) 



1 - iA 2 r 2 



(At s « 1) 



(6) 



which has the same form as Eq. 0, but with important 
differences. First, the memory now rephases only at in- 
teger multiples of r s [32|, [3J, |S3| . Second, the wider is the 
original spectrum n(uS), the larger number of peaks at uik 
with rifc > we can select. For a fixed number of spins 
N' , this allows for narrower A, leading to higher fidelity 
at t s - Alternatively, for the same fidelity (fixed A), this 
results in a larger number of spins N'. 



Our goal is to filter out of the SE with broad spectrum 
n(u)) a subensemble with the spectrum n'{ui) that will 
maximize the resulting fidelity F' while still containing 
many spins N' = fn'(u>)du) ^> 1. Before discussing the 
filtering procedure, let us deduce the optimal spectrum 



III. FILTERING OF THE ENSEMBLE 

We now present the filtering procedure in an ensemble 
of N active dopants, while its optimization is the subject 
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FIG. 1. Level scheme of dopants: \g) and |e) are the ground 
and excited states of the spin-| subspace, and \a) is the aux- 
iliary metastable state. The subensemble with desired spec- 
trum n'{uj) is selected by the il(t) field, which is followed by 
transferring all the remaining dopants to \a). 



of the following Section. 

To select a subensemble of spins with spectrum n'(w), 
we employ another (auxiliary) long-lived state \a) outside 
the spin-i subspace { \g), |e)} of the dopants (see Fig-HJ). 
The preparation of the subensemble proceeds in three 
steps: (i) Starting with all the dopants in the ground 
state \g), apply an external pulse of Rabi frequency Q(t) 
that excites them to state |e). The duration T of the 
pulse should be long enough, in order to select only the 
dopants with transition frequencies ui within a range of 
A ~ 2ir/T around the desired frequency uq, while the 
shape of fl(t) is designed to optimize the resulting fre- 
quency spectrum, (ii) Transfer all the dopants remaining 
in \g) to the auxiliary state \a) by another strong field, 
using an adiabatic sweep across the \g) — > \a) transition 
j4ll 1-421 ] . (iii) Return the dopants selected in step (i) from 
|e) to \g) by, e.g., the adiabatic transfer. 

The chosen subensemble is now ready to use. Its spec- 
trum is given by n'(w) = n{uj)P{uj — luq). where 



P(«-w )= (e\T + e~ l f t > H ^ dt \g) 



(7) 



is the probability to excite the spin with transition fre- 
quency to by the preparation pulse Q(t) used in step 
(i). We consider only the amplitude-modulated field il(t) 
with the fixed carrier frequency ujq', the corresponding 
Hamiltonian is H{t) = ^(w — Ct>o)cr 2 + fl(t)a x with a zx 
the Pauli spin operators. With the new spectrum n'(u>), 
the memory fidelity is given by 



v^ 7 



1 
N 



n(u)P(oj - u^e-^-^^doj. (8) 



Ideally, we would like all the spins at the resonant fre- 
quency w = loo to be selected, P(0) = 1, setting the pulse 
area A = J Vt(t)dt = tt/2 (a 7r-pulse). This still leaves us 
the freedom to choose the shape of f2(£) so as to maximize 
the fidelity in Eq. flgj}. 

As an example, consider a square preparation 
pulse Q(t) = tt/2T (t e [0,T]), for which 
Eq. ([7]) can be be solved exactly, P(w — Wo) = 

^sinc^iV 71 " 2 + (w -cj ) 2 T 2 ). Using Eq. ©, we find 
linear fidelity loss F'(t) sa 1 — 4r/T for short times 
r <C T. With such a "nai've" choice of O(t), the spec- 
tral variance ((u> — ^o) 2 ) docs not converge due to the 



long wings of P(w 
fidelity (Fig. El). 



ujq) (see Fig. [2j, leading to poor 



IV. OPTIMIZING THE PREPARATION 

The ideal rectangular spectrum of subensemble found 
in Sec. lIIBl would require infinitely long preparation time. 
Our goal is therefore to find the optimal preparation 
pulse f2(i) of total duration T which should be shorter 
than the spin relaxation time. To this end, we employ 
the method of Lagrange multipliers to maximize the fi- 
delity F'(t) for a given number of selected spins N' . [Re- 
call that N' determines the transfer time r tr oc (TV') -1 / 2 , 
which in turn should to be smaller than the decay time of 
the (cavity) field.] For convenience, we will actually max- 
imize N'^/F^ using the number spins N' and the pulse 
area A as the constraints. The resulting Euler-Lagrange 
equation reads 



<9[7V , V / P 7 (t)] _ ON' 



dfl(t) 



A, 



+ A 



dA 



9fi(t) d£l{ty 



(9) 



where Ai and A2 are the Lagrange multipliers. We will 
pursue solutions of this equation yielding the optimal 
preparation pulse 51 (t). 



A. Approximate analytic solutions 

Although Eq. (jHl) can be studied numerically, it is in- 
structive to solve an approximate version of this equation 
analytically. 

From Eq. ([?])■ in second order in f2, we have 



P( W ) 



n(t) e -<"-"o)t dt 



(10) 



Assuming that P(uj) is much narrower than the initial 
spectral distribution n(cj), we obtain 



F'(r) 



47T 2 rio 



N 



l2 



n(t + r)n(t)dt 



(ii) 



with TV' sa 2irno L £l 2 (t)dt and no = n{u>o). Note that 
for r > T the fidelity vanishes, i.e., one cannot store QI 
for time r longer than the preparation time T. Equa- 
tion © now reduces to 



Q(t + t) + Q(t -t) = Aifi(i) + A 2 , 



(12) 



where we used the constraint A 
Clearly, near the resonance 



£n(t)dt. 



UJ 



Wo I 



< 



max[Q(t)] 



Eq. (fTO)) is incorrect as the selection probability P(w 
w o) ^ x becomes larger than 1 [see Fig.^b)]. This leads 
to the overestimate of the number of selected spins N', 
making the corresponding constraint inexact. We will sec 
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FIG. 2. (a) Square (blue dashed line) and optimal (black 
solid line) preparation pulses, (b) The corresponding selec- 
tion spectra P(w) calculated numerically. The approximate 
selection spectrum of Eq. (|10[) for the optimal pulse is also 
shown (red doted line); the exact and approximate spectra 
differ significantly around resonance but start to coincide for 
luT/2tt > 1. (c) Second moments uj 2 P(lu) for the correspond- 
ing selection spectra. 




FIG. 3. Fidelity F (r) for subensemble selected by square 
preparation pulse (blue dashed line) and by optimal pulse 
of Eq. (|14p: exact solution (black solid line) and approximate 
analytical solution of Eq. (|15fl (red dotted line) . The duration 
of the pulses is T = 10r tr . Inset: Rabi oscillations of a single 
excitation between the field and the correspondingly selected 
subensemble. 



below, however, that the maximal fidelity obtained un- 
der this approximation is close to the exact, numerically 
calculated fidelity, especially for r <C T [Fig. |3J. This 
is because the loss of fidelity is mainly due the wings of 
the selected spectrum, where the behavior of Eq. (|Tu| is 
correct. We also note that moderate nonuniformity of 
preparation pulse 0(t) for different spins would decrease 
P(u ~ Wo) and thereby the number of selected spins TV', 
but it would affect little the wings of the selected spec- 
trum P(uj) and the resulting fidelity F' . 

We now seek the optimal preparation pulse which will 
maximize the fidelity (a) over a continuous time interval 
t £ [0,Tt r ] (for QI transfer), and (6) at a specific time t s 
(for QI storage). 

(a) To maximize the fidelity for all times r < r tr -C T, 
we notice that Eq. (JT3J) becomes independent of r: 



h(t) = -AiO(t) +A 2 



(13) 



where Ai = — (Ai — 2)/t 2 and A2 = A2/T 2 arc the 
rescaled Lagrange multipliers. The highest fidelity is then 
achieved with the pulse 



il(t) = ft sin(7rf/T), 



(14) 



where Q = f^, so that A = ir/2 and TV' ra 7r 5 n /16T. 
The optimal pulse (|14p and the corresponding (exact and 
approximate) selection spectrum P(u>) and its second 
moment ui 2 P(uj) are shown in Fig. [2] Due to the sup- 
pressed wings of P(w), the spectral variance converges to 
((w — wo) 2 ) ~ 7r 2 /T 2 , while it does not converge for the 
square preparation pulse. The resulting fidelity is shown 



in Fig. |3] and is given by 

(T - t) cos(ttt/T) 



F'(r) 



sin(7rr/T) 



T 



1 



T 2 ■ 

(15) 
Hence, for short times r <C T, the fidelity loss is 
quadratic in t, which should be contrasted with linear 
fidelity loss for the square preparation pulse. 

In the inset of Fig. [3l we show the Rabi oscillations 
of a single excitation between the field and the selected 
subensemble. Using the fidelity of Eq. (|T5l) . we find the 
following solution of Eq. © up to second order in r tI /T: 



a{r) 



1 



5I ) cos(7rr/r t r) 



# (16) 



which yields |a(r tr )| 2 w 1 - 4(r tr /r) 2 and |a(2r tr )| 2 « 
1 + 0(r tr /T) 3 . Remarkably, the probability of excitation 
retrieval into the field mode is higher at the end of the 
second oscillation at t = 2r tr than the end of the first at 
t = T tr- [Note that at the end of the second Rabi cycle 
the combined state of the field and SE has the initial 
phase, |4 f (2r tr )) ~ |1,t/> ).] In general, for this spectrum 
the retrieval infidelity at even revivals is 3rd order in 
Ttr/T, as opposed to 2nd order at odd revivals. 

(b) To maximize the fidelity at a specific storage time 
r s , we assume that the preparation time is a multiple of 
t s , T = mr s (m £ Z). The solution of Eq. (T5J is then 



n{t) = n ^t^T s [t/r s \)s\ 



Wt s \ + 1 



1 



where CIq = -| tan I ^-j- ) , ['"J denotes the integer part 
(floor) of the expression, and £(t) is a temporal profile 
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FIG. 4. Left column: preparation fields Q(t) of Eq. (fl"7j) . op- 
timized for storage time r s , consisting of (a) series of instan- 
taneous pulses £(£) = 8(t — T s /2), (b) series of short pulses 
of width t s /5, and (c) step function constant during each r a 
interval £(£) = l/r s . The preparation time is T = 5r s . Right 
column: the resulting selection spectrum P(td) of Eq. (|TSJ) 
given by the product of comb and envelope. All the spectra 
have the same fidelity F'(t s ) but different number of spins 



within t s of unity area f* £(t)dt = 1. Hence, the prepa- 
ration pulse is an to- fold repetition of £(£), at each step 
< n < to multiplied by a different amplitude, 



£l(t) = ilo £(t—riT s ) sin 7r 



n + 1 

TO+ 1 



fort g [nr 8 ,(n+l)r s ), 

(17) 
so that A = tt/2 and AT' = 7rn (m+l)O§ /J" 3 £ 2 (t)d£. The 
resulting approximate selection spectrum of Eq. (|10l) , 



P( W ) 



n=0 



77 + 1 

771 + 1 



l£MI 2 , (18) 



is a product of two terms: a "comb" term, which is a 
series of optimally shaped peaks spectrally separated by 
2tt/t s ; and an "envelope" term, £(w) = /J" 3 e~ iult (,(t)dt, 
wider than 2tt/t s . Examples of P(w) for various £(£) are 
shown in Fig.|U For a uniform £(i) = l/r s , P(u>) reduces 
to a single peak, since the envelope has a width of 2tt/t s , 
resulting in the minimal number of selected spins N' oc 
1/t s . By contrast, choosing £(£) to be a narrow pulse 
of width St < t 8 yields multiple (~ T s /5t) peaks within 
the envelope, and hence larger number of selected spins, 
N' oc 1/St. A delta function, £(t) = S(t - r/2), will give 
rise to infinitely many peaks and thus an infinite number 
of spins N' — > co. In practice, however, the number of 
selected peaks is limited by the width of the original SE 
spectrum n(u>). The resulting fidelity [Eq. (|lll) ] at time 
r s is now given by 



F'(t s ) 



cos 



T/t s + 1 



(19) 
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FIG. 5. Left column: preparation fields $!(£) optimized for 
(a) long and (b) short preparation time T. In (a) the dashed 
(red) curve is the truncated sinc(Tj) function, and the solid 
(black) curve is the numerically optimized pulse. In (b) the 
dashed (red) curve is the sin(t) function of Eq. (|14[) and the 
solid (black) curve is the numerically optimized pulse. Right 
column: the resulting exact selection spectra P(uj). 



Note that the choice of £(£) does not affect the fidelity 
[unlike the number of selected spins N' oc f Q s £ 2 (t)dt], 
which is to be expected, since different solutions must, 
by definition, yield the same maximal fidelity. 



B. Numeric solutions 

We now solve the Euler-Lagrange equation (J9j> numer- 
ically using the exact selection spectrum of Eq. 0. Wc 
maximize the subensemble fidelity F'(t s ) for two distinct 
situations: 

a. Stable spins - In this scenario, we can use a long 
preparation time T so as to achieve the optimal rect- 
angular spectrum P(w) of the subensemble. Our main 
constraint is the number of selected spins N'. Clearly, 
for infinite T the preparation pulse f2(t) is a sinc(i) func- 
tion; for long but finite T, however, the optimal pulse 
is modified but still resembles the sinc(t), as shown in 
Fig. (Ha). 

b. Short-lived spins - Now the short lifetime of the 
spins severely limits the preparation time T. The result- 
ing spectral width of the selected subensemble A ~ 2tt/T 
will be wide enough to contain many spins. We can there- 
fore disregard the constraint on N', focusing instead on 
achieving the narrowest possible selection spectrum P(w) 
within a given preparation time T. We then find that the 
optimal pulse f2(i) is almost identical to that in Eq. (|14[) . 
see Fig. [5jb). 



C. Fidelity loss during the QI transfer and storage 

The preparation selects a subensemble of N' ^l N 
spins with reduced spectral variance ((w — Wo) ) — 



(^) (N' /no) 2 , but also results in a longer transfer time 
T tr = n/fjy/N'. During the transfer, the fidelity loss, or 
error due to the decay k of the (cavity) field is e tr ~ Kr tr , 
while the error accumulated during the storage time t s 
in the subensemble is e s ~ ((w — aJo) 2 ) T s 2 - Minimizing the 
total error e = e tr + e s with respect to N', we obtain 



\VnaJ 2 4 V" T s 



2\2/5 




? (20) 

Hence, small error e requires large coopcrativity t] 2 uq 3> 
k 2 t s of the field-subensemble coupling. 



V. EXPERIMENTAL CONSIDERATIONS 

To illustrate the results of the foregoing discussion, we 
consider a SE of NV color centers in diamond coupled to 
a superconducting coplanar waveguide resonator [l|-|5J|. 
The ground \g), excited |e) and auxiliary \a) states cor- 
respond, respectively, to the m = 0, m = 1 and m = — 1 
Zeeman sublevels of the ground electronic (spin-triplet) 
state of the NV. Transitions \g) — > |e), \a) have the fre- 
quencies of around 2.88 GHz, and can be selectively ad- 
dressed by the external ct± -polarized microwave fields. In 
addition, a static magnetic field can be used to tune the 
transition \g) — > |e) in and out of resonance with the cav- 
ity mode u>o- The experimental inhomogeneous spectrum 
of the ensemble of N ps 10 12 NV centers [f| has the total 
width of about A/2tt ~ 7 MHz, composed of three par- 
tially overlapping Lorentzians of widths ~ 2.6 MHz split 
by ~ 2.2 MHz due to the hyperfine coupling to the 1=1 
nuclear spin of the 14 N atom. This results in a storage 
fidelity of F(t s ) w 1 — r s /60ns, while with the collective 
SE-cavity coupling strength fjy/N ~ 2ir x 13 MHz Q, the 
excitation transfer time is r tr ~ 40 ns. 

As an example, assume that one can achieve a cav- 
ity quality factor of Q = 10 6 . The photon lifetime 
in the cavity, k~ 1 rj 55 us, is then much longer than 
the transfer time T tr . This allows a preparation that 
reduces the ensemble spectral width A, and the num- 
ber of spins N, by a factor of 5 • 10 3 — almost 4 or- 
ders of magnitude — while still keeping the transfer time 
Tt r ~ 2.8 /is = 0.05k -1 -C k^ 1 well within the cavity 
lifetime. The new subensemble, created by the optimal 
preparation pulse fl(t) of T = 0.7 ms duration, has the 
storage fidelity F'(t s ) « 1 — (r s /0.22ms) 2 . There is neg- 
ligible loss of fidelity during the QI transfer from the 
cavity field to the new ensemble, which can now store QI 
for r s ~ 50 /is with 95% fidelity, compared to 3 ns with 
the original SE and 2.8 us inside the cavity. 

Had we used the square preparation pulse of the same 
duration, the new ensemble would have had a storage 



fidelity of F'(t s ) w 1 — r s /0.17ms, which can store QI for 
r s ~ 8.5 fis with 95% fidelity — more than 5 times worse 
than with the optimal preparation. 

VI. SUMMARY AND DISCUSSION 

Ensembles of long-lived two-level systems — spins — can 
serve as collective quantum memories, but their utility is 
often compromised by the inhomogeneous spectral broad- 
ening, which can hamper both the QI transfer to and 
from the SE and reduce the storage fidelity The fidelity 
is very sensitive to not only the width A but also to 
the profile of the inhomogeneous spectrum of the SE. 
Specifically, for a spectrum with long wings, such as a 
Lorentzian, the loss of fidelity during the storage time 
t s <C A -1 scales linearly with Ar s , while for a spectrum 
with sharp cutoff at A the loss of fidelity is quadratic in 
Ar s . 

In this paper, we have proposed and analyzed a method 
to select a spectrally narrow subensemble of spins, which 
maximizes the fidelity of QI storage and also improves the 
QI transfer from and to an electromagnetic field. In our 
method, judiciously designed pulses of finite duration — 
determined by the spin relaxation time — select the opti- 
mal spectral portion of a large SE, while the remaining 
spins making up the spurious part of the spectrum are 
discarded by transferring them to an auxiliary metastable 
state. Our method is applicable to microwave cavity (cir- 
cuit QED) based hybrid quantum systems involving QI 
processing qubits and SE quantum memories, as well as 
to optical field storage in SEs. In the former case, the 
QI transfer time is limited by the cavity field relaxation 
time, while in the latter case, it is the interaction (tran- 
sit) time of the optical pulse with the SE. Hence, one 
always has to attain large coopcrativity in enscmblc-ficld 
coupling by having many spins, and large optical depth. 

Similar considerations may apply to the schemes in- 
volving noisy processing qubits coupled directly to the 
spin-ensemble memories, such as, e.g., the electron spin 
of a NV center interacting with the surrounding ensem- 
ble of long-lived nuclear spins [44H491 ] . Then in Eq. (|20|) 
one would have to replace k by the qubit decohcrcnce 
rate 7. As opposed to the cavity decay rate k, which is 
difficult to change, 7 of a qubit can be suppressed by dy- 
namical control methods [50H55J , which, for a fixed error 
e, will result in quadratic increase of the memory time 
r s oc e 5 / 2 / 7 2 . 
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